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In this article, we study the global dynamics of a discrete two dimensional competition model. We give sufficient 
conditions on the persistence of one species and the existence of local asymptotically stable interior period-2 
orbit for this system. Moreover, we show that for a certain parameter range, there exists a compact interior 
attract or that attracts all interior points except a Lebesgue measure zero set. This result gives a weaker form of 
coexistence which is referred to as relative permanence. This new concept of coexistence combined with numerical 
simulations strongly suggests that the basin of attraction of the locally asymptotically stable interior period-2 
orbit is an infinite union of connected components. This idea may apply to many other ecological models. Finally, 
we discuss the generic dynamical structure that gives relative permanence. 

Keywords: Basin of Attraction, Period-2 Orbit, Uniformly Persistent, Permanence, Relative Permanence 
AMS Subject Classification: Primary 37B25; 39A11; 54H20; Secondary 92D25 



1. A discrete two species competition model 

Mathematical models can provide important insight into the general conditions that permit the 
coexistence of competing species and the situations that lead to competitive exclusion (Elaydi 
and Yakubu 2002). A model of resource- mediated competition between two competing species 
can be described as follows (Adler 1990; Franke and Yakubu 1991a, 1991b) 

Xn+l = — (1) 

a + Xn + Vn 

Vn+l = yne'-^-(""+^") (2) 

where Xn and i/n denote the population sizes of two competing species x and y at generation n 
respectively; all parameters ri,r2,a are strictly positive. Franke and Yakubu (1991a) established 
the ecological principle of mutual exclusion as a mathematical theorem in a general discrete 
two-species competition system including (l)-(2). They (1991b) also gave an example that such 
exclusion principle fails where two species can coexist through a locally stable period-2 orbit. This 
phenomenon of coexistence has been observed in many other competition models (e.g., Yakubu 
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1995, 1998; Elaydi and Yakubu 2002a, 2002b) including system (l)-(2) with a = 0: 



riXn 



(3) 



pr2-(Xr^+yr^) 



Vn+l Vne 



(4) 



Notice that the equation (3) is the non-overlapping lottery model (Chesson 1981) with singularity 
at the origin. Every initial condition (xo,0) with xq > maps to (ri,0). The lottery model 
emphasizes the role of chance. It assumes that resources are captured at random by recruits from 
a larger pool of potential colonists (Sale 1978; chapter 18, Chain et al 2011). When yn = 0, (3) 
can be a reasonably good approximation for plant species where a single individual can sometimes 
grow very big in the absence of competition from others or for a territorial marine species, such 
as coral reef fish, where a single individual puts out huge number of larvae (communications with 
P. Chesson; also see Sale 1978). Chesson and Warner (1981) used such non-overlapping lottery 
models to study competition of species in a temporally varying environment. In this paper, we 
focus on the dynamics of (3)- (4). The system (3)- (4) may be an appropriate model for resource 
competition between a territorial species x and a non-territorial species y. 

A recent study by Kang (submitted to JDEA) shows that (3)- (4) is persistent with respect to 
the total population of two species, i.e., all initial conditions in \ {(0,0)} are attracted to a 
compact set which is bounded away from the origin. The results obtained in Kang (preprint 2010) 
allow us to explore the structure of the basin of attraction of the asymptotically stable period-2 
orbit of the system (3)- (4) lying in the interior of the quadrant. In this article, we study the global 
dynamics of (3)- (4). The objectives of our study are two- fold: 

(1) Mathematically, it is interesting to study the global dynamics of (3)-(4) since it has sin- 
gularity at the origin. Thus, the first objective of our study is to give sufficient conditions 
for competitive exclusion and coexistence of (3)-(4). 

(2) Biologically, it is very important to classify and give sufficient conditions for the coexis- 
tence of species in ecological models. Among many forms of coexistence, permanence is 
the strongest concept since it requires all strictly positive initial conditions converge to 
the bounded interior attractor. Although permanence fails for (3)- (4), we establish the 
weaker notion relative permanence: almost all (relative to Lebesgue measure) strictly pos- 
itive initial conditions converge to the bounded interior attractor. Numerical simulations 
of other ecological models (e.g., Franke and Yakubu 1991; Kon 2006; Cushing, Henson and 
Blackburn 2007; Kuang and Chesson 2008) suggest the possibility that relative perma- 
nence may apply where permanence fails. Our second objective of this article is to draw 
attentions on the concept of relative permanence. Our study could potentially provide 
insight on weaker forms of coexistence for general ecological models and open problems 
on the basins of attractions of stable cycles for a discrete competition model studied by 
Elaydi and Yakubu (2002a). 

Simple analysis combined with numerical simulations suggest the following interesting dynamics 
of the system (3)-(4) 

(1) There is no interior fixed point. The eigenvalue governing the local transverse stability of 
the boundary equilibrium on the x-axes (i.e., y = 0) is given by e^^-^i (^or ^ on the y- 
axes). If this eigenvalue is less than 1, then we say that the equilibrium point on the x-axes 
(or ?/-axes) is transversally stable, otherwise, it is transversally unstable. Thus, if ri > r2, 
then the boundary equilibrium ^* = (^1,0) is transversally stable and 77* = (0,r2) is 
transversally unstable; while ri < r2, ^* = (ri, 0) is transversally unstable and 77* = (0, r2) 
is transversally stable. 

(2) For a certain range of ri and r2 values, there exists an asymptotically stable periodic-2 
orbit in the interior of the quadrant which attracts almost every interior point in R^. For 
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example, when ri — 2^r2 — 2.2, the periodic-2 orbit is given by 

{x\,y\) = (0.1536,2.9629) and {x\,y\) = (0.0986,1.1849) 

and the eigenvalues of the product of the Jacobian matrices along the orbit are 0.91 and 
0.26. 

(3) There exists a heteroclinic orbit connecting ^* to 77* (see Figure 1); 



A heteroclinic orbit for a discrete 2-species competition model 



r^=2, x^=2.2, X|j=2, y^jsCOOl, 
Iterations N=100 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

Species X 



Figure 1. A heteroclinic orbit of the system (3)-(4) when ri = 2, r2 = 2.2, = 2, yo = 0.001. 

(4) The basin of attraction of the interior periodic-2 orbit P2 consists of all interior points of 
except all the pre-images of the heteroclinic curve C (where C is the closure of the 
union of all heteroclinic orbits, see Figure 2). 

Basins of attraction of P' 




Figure 2. The basin of attraction of the interior period-2 orbit is the open quadrant minus the pre-images of the heteroclinic 
curve C . The latter partition the quadrant into components which are colored according to which of the two periodic points 
attract points in the component under the second iterate of the map. Given a point in one of the regions, there is a large 
number A^, such that the point will be very close to {x\^y\) at the iteration t and will be very close to (x2, ^2) iteration 
f + 1 for ah t > N. 



Moreover, further analysis and numerical simulations suggest that if the system (3)- (4) satisfies 
the following conditions C1-C3, then it has the same global dynamics as the case ri = 2, r2 = 2.2: 



• CI: The values of ri,r2 satisfy 

2 < r2 < 2.52, r2 > ri > 1, and e^^^"^"^''" 



> 1. 



• C2: There is a boundary period-2 orbit My = {771,772} = {{0,yi), (0,^2)} where > 1. 

• C3: There is a heteroclinic orbit connecting ^* to 77* (see Figure 1). 
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Condition CI implies that the equihbria ^* and 77* of the system (3)- (4) are saddle nodes, where 
^* is transversally unstable and 77* is transversally stable. Moreover, species y can invade species 
X. Condition 2 < r2 < 2.52 combined with Condition C2 indicates that species x can invade 
species y on its periodic-2 orbit {(0,?/i), (0,^2)}- Figure 3 describes the schematic scheme of the 
global dynamics of the system (3)-(4) when it satisfies Condition C1-C3. 




Figure 3. Schematic features of the system (3)-(4) when ri = 2.,r2 = 2.2. 



The structure of the rest of the article is as follows: In Section 2, we give the basic notations 
and preliminary results that will be used in proving our main theorems. In Section 3, we obtain 
sufficient conditions on the persistence of one species and the extinction of the other species by 
using Lyapunov functions (Theorem 3.1). In Section 4 we first give a sufficient condition on the 
existence of locally asymptotically stable interior period-2 orbit for the system (3)- (4) (Theorem 
4.1); then we show that for a certain parameter range, the system (3)- (4) is relative permanent, i.e., 
it has a compact interior attractor that attracts almost points in (Theorem 4.3) by applying 
theorems from persistent theory. In the last section 5 we discuss the fact that the global dynamics 
of the system (3)- (4) are generic rather than rare. Similar dynamic behaviors of (3)- (4) have been 
observed in many biological models. Studying sufficient conditions for the relative permanence of 
the generalization of such biological models can be our future direction. 



2. Notations and preliminarily results 

Notice that the system (3)-(4) has singularity at the origin (0, 0), thus its state space is defined as 
X = {(x, y) <x + y < oc}. Let H denote the map defined by (3)-(4). Then H : X ^ X 

is a discrete semi-dynamical system where H^{^q) = Co = (^o, Vo) and H^{^o) = = (^n, yn)iT^ ^ 
Z+. Here, we give some definitions that will be used in the rest of the article. 

Definition 2.1. [Pre-images of a Point] For a given point ^ we say ^ G X is a rank-fc 
pre-image of ^0 if H^{i) = ^0- The collection of rank-fc {k > 1) pre-images of is defined as 

and the collection of all pre-images of (including /c = 0) is defined as 

EF^,^l\jH-H^o)]\J{^o}. 
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Definition 2.2. [Invariant Set] We say M C X is an invariant set of H if H{M) — M. 

Definition 2.3. [Pre-images of an Invariant Set] Let M be an invariant set for the system 
(3)-(4), then H^{M) = H{M) = M. The collection of rank-/c pre-images of M (A: > 1) is defined 

as 

H-\M) = [J{^eX\M: H\0 = Co}; 
and the collection of all pre-images of M (including A: = 0) is defined as 



EFm= (uk>iH-''{M))[jM 



U U {^eX\M:H''iO=^o} 



k>i UoeM 



[JM. 



Note: If M is an invariant set of i7, then H ^{M) should not contain points in M for all /c > 1. 

Definition 2.4. [Uniform Weak Repeller] Let X be a positively invariant subset of X. We call 
the compact invariant set C is a uniformly weak repeller with respect to X if there exists some 
e > such that 



limsupd(i?^(0,C) > e for any ^ G X \ C. 



Definition 2.5. [Uniform Weak p-Persistence] Let X be a positively invariant subset of X. The 
semi-flow H is called uniformly weakly p-persistent in X if there exists some e > such that 

limsup/}(i7'^(^)) > e for any ^ E X 

n— T^oo 

where p : X ^ is a persistence function (e.g., p(x^y) = x can be a persistence function if 
we want to study whether species x is uniformly weakly persistent or not). We say species x is 
uniformly weakly persistent in X if there exists a e > such that 

limsupxri > e for any ^ G X. 

Definition 2.6. [Uniform Persistence] Let X be a positively invariant subset of X. We say 
species x is uniformly persistent in X if there exists some e > such that 

liminf Xn > e for any ^ G X. 

Definition 2.7. [Permanence] Let X be a positively invariant subset of X. We say the system 
H is permanent in X if there exists some e > such that 

liminf min{x^, y^} > e for any ^ G X. 

Definition 2.8. [Relative Permanence] We say the system H is relative permanent in X if there 
exists some e > such that liminf^^oo i^ii^l^n? I/n} > ^ foi" almost all initial condition taken in 
X (i.e., all initial conditions in X except a Lebesgue measure zero set). 

Lemma 2.9. [Compact Positively Invariant Set] Assume that ri ^ r2, then for any 

< e < min{ri,r2,e2^^-^-^^'"',rie^^-^^} = r^. 
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the compact region defined by 

De = {{x,y) e X : e < x + y < max{ri, e''""^}} 

is positively invariant and attracts all points in X . 

Lemma 2.10. [Pre-images of Invariant Smooth Curve] Assume that 

r2 > ri > 1 and e^^^-^-^^'"' > 1. 

Let C he an invariant smooth curve of the system (3) -(4) and M he any compact suhset of X , 
then 1712 {EFc H M) = where m2 is a Lehesgue measure in 

Remark: Lemma 2.9 and 2.10 are a direct corollary from Theorem 2.2 and 3.3 in Kang (sub- 
mitted to JDEA). 

3. Sufficient conditions for persistence 

In this section we investigate sufficient conditions for the extinction of one species and the persis- 
tence of the other species in system (3)-(4). Let be the set defined in Lemma 2.9 and denote 
D(: as the interior of D^. We can obtain sufficient conditions for the extinction of one species by 
using a Lyapunov function V : ^ where V{x, y) — x^y^ and c, d are some constants. 
In addition, we give a sufficient condition on the persistence of species y by applying Theorem 
2.2 and its corollary of Hutson (Hutson 1984) through defining an average Lyapunov function 
P{x^y) — y in the compact positively invariant region D^. Now we are going to give detailed 
proof of the following theorem: 

Theorem 3.1 . [Persistence of One Species] 

(1) If ri > r2 > 0, then the system (3) -(4) has glohal stahility at (ri,0); i.e., for any initial 
condition = (xo^yo) G {(x^y) G X : xo > 0}^ we have 

lim i/"(Co) = lim H'^ix^y) = (ri,0). 

(2) If < ri < r2, then the species y is uniformly persistent in X, i.e., there exists a positive 
numher 5 > such that for any initial condition — (^o^^/o) ^ {{^^u) ^ ^ - U > 0}; 
have 

lim inf yn ^ ^ 

where {xn^yn) — H^i^o)- Moreover, if ^^"^^ ^ > 1, then the species x goes to extinct 



ri 

for any = (^o, Vo) ^ {{x, y) e X :y >0}, i.e.. 



lim Xn = 0. 



Proof. According to Lemma 2.9, any point in X is attracted to the compact positively invariant 
set for any e G (0,r^]. Therefore, we can restrict the dynamics of (3)-(4) to D^. 
If ri > r2 > 0, define V{x, y) = x~'^^y, then 

V{Hix,y)) , 

..\ri ^r2-x-y 



Vix,y) 
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Let f{u) = r^^^u'^^e'^^~^. Since f^{u) = r^^^u'^^~^{ri—u)e'^'^~^, we can conclude that the maximum 
value of f{u) achieves at 7i = 1, i.e., 

max {f{u)} = /(ri) = (^!l)^ig^2-ri ^ ^ where K > max{ri, r2, e^^""*^, 1}. 

e<u<K r\ 

According to Lemma 2.9, we know that is positively invariant and attracts all points in X. 
Therefore, any point in the region {(xq, yo) ^ : xq > 0} has the following two situations 

(1) If yo = 0, then V (H^(xo, yo)) = V ((x,, y^)) = or 

(2) If xo, yo > 0, then 

y(g(xo,^o)) ^ niax K-(x + ^re-— n 

= max {rr^^ix^^e^^-^l = (!l)^ie^^-^^ < 1. 

e<u<K ^ n 

Thus, 



V (i/-(xo, i/o)) _ V (i/-(xo, yo)) (i?(xo, yo)) 



V(xo,yo) 



y^-i(xo,yo) 



^(^o,yo) 



as n ^ oc. 



Therefore, the positively invariant property of implies that 



Therefore, 



This indicates that 



lim x^'^'yn 



lim yn — ^ and liminf > 6. 



lim Xn+i — lini = lim — 



Therefore, if t\ > T2 > 0, then the system (3)-(4) has global stability at ^* = (^i,0). The first 
part of Theorem 3.1 holds. 

If T2> r\> 0, then the omega limit set of Sx = {(x, 0) : x > 0} is i.e., ijj{Sx) = {C*}- The 
external Lyapunov exponent of Sx is e^^"^^ > 1, therefore, it is transversal unstable. According 
to Lemma 2.9, for any < e < r^, attracts all points in X . Thus, the uniform persistence 
of species y follows from Theorem 2.2 and its corollary of Hut son (Hutson 1984) by defining a 
Lyapunov function P(x^y) = y on the compact positively invariant region D^, i.e., there exists a 
positive number 6 > such that for any > 0, we have 

lim inf y^ > 6. 

If, in addition, ri < e^^2-i-e'^2 1^ ^j^gj^ ^^^^^ define a Lyapunov function as V{x, y) = xy~^^ then 
we have 



ri 



V{H{x,y)) ^ 

V{x,y) (x + y)e^2-(x+2/)" 
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Now choose e = min{ri, r2, e^^^"-*^"^"^^ \^ie^^~^^}, then e = ri since r2 < ri < e^^^-i-e'^s 
Therefore, any point (xq, ^o) ^ satisfies ri < x + y < e^^"-^ and will stay in D^^ for all future 
time. This implies that for any point (xq, yo) in Dr^ with yo > 0, we have 

V{H{xo,yo)) ^ ri ri 



V^(xo, yo) {x,y)eD^, {x + y)e^2-(^+^) min(^^^)^^^^ {(x + y)e^^-i^+y)} 

ri ri 



l-e-2 



T<1 



Hence, lim^^oo^n — 0. Now if (xo,^o) ^ \ ^ri? then according to Lemma 2.9, (xo,^o) will 
either enter D^^ in some finite time or converge to (ri,0). Now we consider the following two 
cases for any initial condition (xq, yo) ^ \ with yo > 0: 

(1) If xo = 0, then x^ = for all positive integer n; 

(2) If Xo > 0, then (x,^) will not converge to (ri,0) since the equilibrium point (ri,0) is a 
saddle and transversal unstable when r2 > ri, therefore, (x, y) will enter Dr^ in some finite 
time. 

Thus, the condition ri < r2 and ri < e^^2-i-e^2 i g^Q^^Q^j^^ggg ^j^^^^ 



lim Xn = 0. 



Therefore, the second part of Theorem 3.1 holds. □ 

Remark: The first part of Theorem 3.1 can be considered as a special case of rational growth rate 
dominating exponential (Franke and Yakubu 1991) which states that if species x with rational 
growth rate can invade species y with exponential growth rate at species y's fixed point, i.e., 
(0,r2) is transversal unstable, then the exponential species goes extinct irrespective of the initial 
population sizes. The second part of Theorem 3.1 shows that the exponential species can persistent 
whenever (ri,0) is transversal unstable (i.e., r2 > ri). However, the rational species may not go 
extinct unless ri < e^^2-i-e"2 ^ ^^^^^ simulations (e.g.. Figure 2) suggest that two species of 
the system (3)- (4) may coexist for almost every initial conditions in X under certain conditions. 
This point will be illustrated with greater details in the next section. 



4. Coexistence of two species 

In this section, we give sufficient conditions for the existence of the interior period-2 orbits and 
its local stability for the system (3)- (4) as the following theorem states: 

Theorem 4.1 . [Sufficient conditions on the existence of interior period-2 orbits] Ifr2 > 2, then 
the Richer map yn+i = yn^^'^~^'' has period two orbits {^1,^2} where < yi < r2 < y2 cind 
yi + y2 — 2r2. The system (3)-(4) has an interior period-2 orbit P2 = {{x\^y\)^ (^2' ^2)} '^here 

i ^ gi(gie"2-^i-g2) i ^ si(s2-ri) 

4 = ^, y2-yW'' (5) 

si x\ + yl = r2 - \/r| - rf, 82 x\ + y| = r2 + ^Jr\ - r\ 
if one of the follows holds 

(1) 516^2-^1 > 82, or 

(2) r2 - y^rl - rj > yi, or 

(3) 2 < n < r2 < 2.5 and ri > r2 - or 
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(4) 2.085 < n < r2 < 2.5 

In particular, (4) implies (3); (3) implies (2) and (2) implies (1). Moreover, if ri — 2 and 
5 — r2 — ri — r2 — 2 is small enough, then P2 is locally asymptotically stable. 

Proof . If r2 > ri > 0, then 



r2 - ri - y r| - = \/r2 ri (\/r2 ri - \/r2 + ri) < 0, 
thus we have the fohowing inequahties: 



51 = r2 - y r| - < ri < r2 < S2 = ^2 + Y ^2 ~ ^i- 

Therefore, from (5), we find that sie^^-^i — 52 > is a sufficient condition for the existence of 

Notice that the Ricker map y^+i — yn^^'^~^"' goes through period-doubling two bifurcation at 
r = 2, thus if r2 > 2, the Ricker map has a period-2 orbit {^1,^2} where 

< ^1 < r2 < ^2, y2 = yie''^~^' and ^1 + ^2 = 2r2. 
Since 5i + 52 = 2r2, then from the graphic representation (see Figure 4), we can see that 



sie 



r2-si 



52 > whenever yi < si < r2. 



Therefore, the condition 5i = r2 — y^r^ — > yi is a sufficient condition for sie^^ _ ^ 
therefore, it is a sufficient condition for the existence of 




Figure 4. The location between Si, i = 1, 2 and sie^2 si , xhe solid line is f{y) = ye^^ ^; the dashed line is f{y) = y; the 
dot line is f{y) = 2r2 — y. 



Let a = Y^rl — r^, then we have the following equivalent relationships: 



sie 



r2-si 



> S2 (r2-a)e'' > r2 + a 



2a(e« + 1) 2a 



e« - 1 



e« - 1 



r2 - a > 



2a 



e« - 1 



(6) 



Thus we find that r2 — a > implies s\e^'^ — S2> 0. 
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If 2 < n < r2 < 2.5, then 

10 



<a = Jrl-rl = - rWr2 + n < \/(L5x/(2.5 + 2.5) 



2 

Notice that h{a) — is a decreasing convex function with respect to a, thus 

2-4^ 



/i(a) < Ha) = 2 J=^a 



2 



o yio 

^ ^/lO 



where A:(a) is a straight hne going through (0,2) and (^^^^ , h{^^)^ . Since — ^^7= — ^ > 0.74, 
therefore, 



2 



^ -, In A/ID 

2 - 0.74a > 2 > for ah < a < 

- v20 - ga _ 1 ' 2 

2 

Hence, from (6), we can conclude that r2 — a > 2 — 0.74a imphes r2 — a > ^Izj, therefore it imphes 
5ie^2-5i _ ^2 > 0. Notice the fohowing equivalent relationships. 



^ o Ci ^ ^ '^2 2 ^ ^ 2 2 ^ 2 2 ^ ^ ^ ( 0^.26 )^ (n\ 

r2-a>2-0.74a^a< — ^r,-r,<{—) ^nyr^-^;^^, (7) 

therefore, we can conclude that ri > r2 ^^4^ implies r2 — a > 2 — 0.74a, therefore, it implies 

r2 — a > therefore, it implies sie^^-^i _ > 0, therefore, it implies the existence of P^- 

Since 2 < ri < r2 < 2.5, then r2 + ri < 2r2 < 5, thus 

{ r2-2 \2 / r2-2 \2 / r2-2 \2 

> ^2 - ^ ri > r2 - ^ n > r2 - ^ ^'^^ ^ 



5 2r2 r2 + ri 

/ 7-2-2 N 2 / r2-2 N2 

Therefore, ri > r2 — ^-^^r — implies ri > r2 — ^^rq?^, therefore, it implies the existence of 
Notice that 

_ ^ Q-26 ^ ^ -2.958579881(r2 - 2.169000001)^ + 2.08450001 < 2.085, 
5 

hence we can conclude that 



2.085 < ri < r2 < 2.5 



( ^2-2 \2 

implies ri > r2 r^'^n ' therefore, it implies the existence of 

So far, we have shown the first part of Theorem 4.1. Now we are going to see that local stability 
of Let ri = 2 and 5 = r2 — ri = r2 — 2, then we have 



si = r2- Jr^-rl = 2 + S- ^/S(4 + S) and 52 = 2 + 5 + ^^(4 + 5). 
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Thus if 6 is small enough, then 

sie^'-'' - 52 = (10/3)5^/2 - (10/3)52 + (41/12)5^/2 ^ Q^gs^ ^ q 

Therefore, from the proof for the first part of Theorem 4.1, we can conclude that the system 
(3)- (4) has an interior period-2 orbit P2 when ri = 2 and 5 = r2 — ri = r2 — 2 is small enough. 
The local stability of P2 is determined by the eigenvalues of the product of the Jacobian matrices 
along the periodic- 2 orbit which can be represented as follows: 

J\p^-\ . I . (8) 

If 5 is small enough, then the trace and determinant of (8) can be approximated by 

det{J) + 1 = 2-^ + ^ + 0{5^) and trace{J) ^2-^ + ^ + 0{5^). 
By the Jury test (p57 in Edelstein-Keshet 2005), we see that P2 is locally asymptotically stable if 

2 > 1 + det{J) = 2 - ^ + ^ + 0(^3) > \trace{J)\ = |2 - ^ + ^ + 0{5')\ (9) 

which is true when 6 is small enough. 

Therefore, the statement of Theorem 4.1 holds. □ 

Remark: Theorem 4.1 provides a sufficient condition on the existence of the interior period 
two orbit and their stability. Numerical simulations suggest that the system (3)- (4) has an interior 
period two orbit whenever 2 < ri < r2 < 2.5. In the case that ri = 2 and r2 = 2 + 5, the interior 
period two orbit P2 is locally asymptotically stable whenever S < 0.95 (see Figure 5-6). 




Figure 5. Interior period two orbit P2 of the system (3)-(4) when ri = 2,r2 = 2 -\- 6 and 6 is varying from to 2. The solid 
line is (x2,y2) and the dashed line is (x\,yl). 



Lemma 4.2. [Pre-images of Heteroclinic orbit] Assume that the system (3)-(4) satisfies Condi- 
tion CI and C3, then there exists a smooth invariant curve C that connecting ^* to 77*. Denote 
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5 

Figure 6. The stability of the interior period-2 orbit P2 of the system (3)-(4) when ri = 2,r2 = 2 -\- S and S is varying from 
to 1. The solid line is det{J) + 1; the dashed-dot line is constant 2 and the dot line is \trace{J)\. This figure indicates that 
P2 is locally asymptotically stable when n = 2 and 2<r2 = 2 + (5< 2.95. 

EFc as the collection of all ranks pre-images of C, then m2(EFc H M) = 0^ where M is any 
compact subset of X and m2 is a Lebesgue measure in R^. 

Proof. First we show that C is a smooth curve connecting ^* = (ri, 0) to 77* = (0, r2). Let 
be the local unstable manifold of ^* and cji{r]*) be the local stable manifold of 77*. Since the map 
H is smooth, then according to stable manifold theorem (Theorem D.l in Appendix, Elaydi 2005), 
we can conclude that both cj'i'iC) '^f (^*) smooth curves. Since (3)- (4) satisfies Condition 
C3, then there exists some positive integer k such that (cl;J^(^*)) is smoothly connected with 
uji{r]''). Thus C is a smooth invariant curve with <^*,77* as its two end points. Then according to 
Lemma 2.10, the statement holds □ 

4-1 Persistence of species x in new space X 
Let 

= {(x, 0) G D,}, Sy = {(0, y) G D,}, S = S,[JSy = {(x, y) e D, : xy = 0}. 

Let C be the closure of all heteroclinic orbits connecting from ^* to 77* and denote EFc as the 
collection of all rank pre-images of C. Then we can conclude that E = EFc H is compact and 
forward invariant. Define the following new spaces 

s = s\Je,x = x\s. 

Then both S and X is positively invariant. In addition, S is compact since both S and E are 
compact. 

The rest of this section, we assume that the system (3)-(4) satisfies Condition C1-C3. We will 
prove the following theorem: 

Theorem 4.3 . [Relative Permanence] Assume that the system (3)-(4) satisfies C1-C3. Denote 
EFc as the collection of all pre-images of the heteroclinic curve C . Then there there exists a 
compact interior attractor in R\ that attracts all points in the interior of X except points in 
EFc' In particular, the interior attractor attracts almost every point with respect to a Lebesgue 
measure in of any compact subset M in the interior of X , i.e., 1712 {EFc fl ^) — 'where 7712 
is a Lebesgue measure in R^. 

Proof. We use the following three main steps to prove the statement. We provide the detailed 
proof of the first two steps in the Appendix and the remaining proof here. 
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(1) C is a uniform weak repeller with respect to X, i.e., there exists some 6 > 0, for any 
^ G X, we have 

hmsupd(i/^(0,C) > b. 

The detailed proof of this part has been shown in the Appendix (Lemma A.l). This imphes 
that any point in X is going to be away from C in some distance in some future time even 
if the point is very close to C, . 

(2) Species x is uniformly weakly persistent in X, i.e., there exists some 5 > 0, such that for 
any initial condition — (^o? Ho) ^ X, the system has 

lim sup Xn > S. 

The detailed proof of this part has been shown in the Appendix (Lemma A. 2). This implies 
that any point in X is going to be away from S in some distance in some future time even 
if the point is very close to S. 

(3) Species x is uniformly persistent in X, i.e., there exists some e > 0, such that for any 
initial condition ^ = (x^y) G X, the system has 

liminf > e. 

Now we will show that the last step. Define a continuous and not identically zero persistent 
function p(^) = (i(^, S) where ^ G X. Then by the definition of the persistent function p, we have 

5 = {^eX:p(if*(0) = 0,Vt>0} 

is nonempty, closed and positively invariant. In addition, the system H satisfies the following two 
conditions: 

(1) There exists no bounded total trajectory (j) such that p{(j){0)) = 0,p((/)(— r)) = and 
p((j)(t)) > for some positive integer r, t. 

(2) H has a compact attractor which attracts all points in X. 

From Lemma A. 2, we know that species x is uniformly weakly persistent in X. Thus by applying 
Theorem 5.2 Smith & Thieme 2010, we can conclude that species x is uniformly persistent in X, 
i.e., there exists some e > 0, such that for any initial condition ^ = (x, y) G X, the system has 

liminf > e. 

Notice that species y is uniformly persistent in X whenever r2 > ri according to Theorem 3.1. 
Thus, species y is also uniformly persistent in X since X is a positively invariant subset of X. 
Therefore, based on the argument above, we can conclude that there exists some positive constant 
e > 0, such that for any initial condition taken in X, we have 

liminf minjxri, ^n} > ^• 

Hence there there exists a compact interior attractor that attracts all points in the interior of X 
except points in EFq. 

Let M be any compact subset of the interior of X. Then any initial condition taken in M 
will enter D^: in some future time through the following two cases 

(1) ^0 ^ EFc which will enter C in some finite time. 
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(2) ^0 ^ which is attracted to the interior compact attractor. 

Since there are only Lebesgue measure zero of points in M that belong to EFq^ therefore, ac- 
cording to Lemma 2.10, we can conclude that m2 {EFcP\M) = for any compact subset M in 
the interior of X. □ 



5. Discussion and future work 



In this article, we study the global dynamics of the system (3)- (4). We give sufficient conditions for 
the uniform persistence of one species and the existence of locally asymptotically stable interior 
period-2 orbits for this system. We also show that for a certain parameter range, the system 
(3)-(4) is relative permanent^ i.e., there exists a compact interior attractor that attracts almost 
all points in X. Numerical simulations strongly suggest that this compact interior attractor is the 
locally asymptotically stable interior period-2 orbit P2 and its basin of attractions consists of a 
infinite union of connected regions that are separated by all pre-images of the heteroclinic curve 
C (see Figure 2). 

The results that we obtained in Theorem 3.1 are a special case for model (l)-(2) when a = 0. 
Our Theorem 4.1 can be extended to the general model (l)-(2) when a > 0. If r2 > 2 and 
r\ > 2ar2, the explicit expressions of the interior periodic-2 orbit {{x\, yl), {x^, y^} of the system 
(l)-(2) can be found as 

i _ (a + gi) (^2 - sie"^-^^) ■ _ risi - S2 (a + ^i) 
ri - (a + 5i)e^2-5i ' ~ _ + si) e''^-^^ 

i _ (a + S2) (si - 826''^-^^) ■ _ riS2 - si (a + S2) 
~ ri - (a + 52)e^2-52 ' ^'^ ~ ^ - {a + S2) e''^-^^' 

where 



51 = x\+y{=r2 + ^ (r2 + a)2 - r^, 52 = 4 + I/2 = ^2 - y (^2 + a)2 - r^. 

This interior periodic-2 orbit can have local stability for a certain range of parameters' values. 
For instance, if 

n = 2.1, a = 0.1 and r2 = 2.5, 
then the system (l)-(2) has locally stable interior periodic-2 orbit 

{xly\) = (0.17,0.80) and = (0.33,3.70) 



along which the eigenvalues of the product of the Jacobian matrices are 0.11 and -0.24. Moreover, 
numerical simulation suggests follows 

(1) There exits a heteroclinic orbit C connecting ^* to 77* (see Figure 7); 

(2) The basin of attraction of the interior periodic-2 orbit P2 is all points in the interior of 

except a Lebesgue measure zero set in which is a collection of all pre-images of the 
heteroclinic curve C (see Figure 8). 

However, more mathematical techniques need to be developed in order to obtain results similar 
to those in Lemma 2.10 for the system (l)-(2) when a > 0. This is an area for future study. 

Our results may apply to a two species discrete-time Lotka-Volterra competition model with 
stocking where both species are governed by Ricker's model and one species is being stocked at 
the constant per capita stocking rate si per generation (lO)-(ll) (Elaydi and Yakubu 2002a&; 
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The heteroclinic orbit connecting from (2,0) to (0,2.5) 



r=2.1, a=0.1, r =2.5, x =2, y =0.00001, 60 iterations 

1 ' '2 ' ' 'O ' 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

Species x 



Figure 7. The heteroclinic orbit of the system (l)-(2) when ri = 2.1, a = 0.1, r2 = 2.5. 




Population of species x 



Figure 8. The basin of attraction of the interior period-2 orbit of the system (l)-(2) when ri = 2.1, a = 0.1, r2 = 2.5 is the 
open quadrant minus all pre-images of the heteroclinic curve C. The latter partition the quadrant into components which 
are colored according to which of the two periodic points attract points in the component under the second iterate of the 
map. Given a point in one of the regions, there is a large number A^, such that the point will be very close to (x\,y\) at the 
iteration t and will be very close to (x2,y2) at the iteration t + 1 for all t > N. 



2002b). We may infer from simulations (see Figure 9) that the basin of attraction of the 2-cycle 
is the infinite union of connected regions that are separated by ah pre-images of the heterochnic 
curve C when si — .5,gi — 1.5,^2 — 2.2, 2?i —p2 — ^- 



Vn+l Vne 



Si + e 

q2-p2(Xr,+yn,) 



ql-pl(Xr^-\-yn) 



(10) 

(11) 



The system (l)-(2) is not the only competition model that has a local stable interior period- 
2 orbit that attracts all points of except all pre-images of the heteroclinic curve C that is 
connecting two nontrivial boundary equilibria. In general, if a discrete two-species competition 
model satisfies the following conditions (see Figure 10 for a schematic presentation), numerical 
simulations suggest that it can have an interior attractor that attracts all points of except all 
pre-images of the heteroclinic orbit that is connecting two nontrivial boundary equilibria. It will 
be our future work to develop more powerful analytic tools to rigorously prove this. 

• Gl:The system has only two nontrivial boundary equilibria (x*, 0) and (0, y*). Moreover, species 
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A phase plane for a competition model 

3r 




Species x 

Figure 9. A single forward orbit of the system (lO)-(ll) starting near the fixed point on x-axis. 

y is persistent. 

• G2:The omega limit set of ?/-axis is a unique attracting period-2 orbit My = {^1,^2} = 
{(0,^1), (0,^2)} on y-axis, which attracts ah points in y-axis except a Lebesgue measure zero 
set. In addition, the external Lyapunov exponent of My is greater than 1, i.e., species x can 
invade species y on My. 

• G3:There is a heteroclinic orbit connecting the boundary equilibrium (x*,0) to (0,7/*). 



y 




X 



Figure 10. The general structure of dynamics that the basin of attraction of interior attractors (e.g., the interior periodic-2 
orbit) is all points in except the collection of all pre-images of the heteroclinic curve C. 
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Appendix A. Important Lemmas 

Lemma A.l. [Uniform Weak Repeller] If the system (3)-(4) satisfied satisfies Condition C1-C3. 
Then there exists some h > Q, such that 

limsupd(i?^(0,C) >h for any i^X. 

Proof . The condition 2 < r2 < 2.52 indicates that species y has a unique attracting period-two 
orbit 

{ni^m} = {(0,yi),(0,^2)} 

in its single state and the condition ri < r2 imphes that the boundary equihbrium 77* = (0, is 
a saddle. By Hartman-Grobman-Cushing theorem (Elaydi 2005), there exists some neighborhood 
of such that any point 77 G Ue^{r]'') nX will exit from this neighborhood in some finite 
time. If we choose e small enough, then 77 is attracted to a compact neighborhood 

B^lMM^j^i^eX :d{^,My)<S} 

where My = {^1,^2} in some finite time. Similarly, the condition < ri < r2 implies that the 
boundary equilibrium ^* is also a saddle, by Hartman-Grobman-Cushing theorem, there exists 
some neighborhood f/eslC*) of C*5 such that any point ^ G Ue^i^"") H X will exit from this neigh- 
borhood in some finite time. 

Choose e = min{ei, 62}. Let K — C \ ((7e(^*) IJ Ue{r]'^))^ then K is a compact subset of C. Since 
C is the closure of the family of heteroclinic orbits connecting ^* to 77*, then any point ^ ^ K will 
reach Ue{r]'') in some finite time me{^). Moreover, there exists a neighborhood of ^, denoted by 
Us^{^) will contain in Ue{r]'') in time 7ne(^), i.e.. 

Then we can see that 

Kc[j Us.iO- 
Since K is compact, it has a finite open cover, i.e., 

fh 

i=l 

Choose S = min{e^mini<i<fn{S^.}}. Then any point ^ G X with d(^^K) < 5, then there exists 
some 771^ = rriei^i), 1 < i < fh, such that H^^{^) G ^7e(^*)- 

Now assume that the statement of Lemma A.l is not true. Then for any k large enough, there 
exists some ^/e G X and a positive integer Uk such that 

< p for any n > n^. (Al) 

Choose k > ^. Then d {H'^'' {^k) 1 C) < S. We show the contradiction in the following three situa- 
tions: 
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(1) If H^^{^^k) ^ Ue{r]*), then by Hartman-Grobman-Cushing theorem, H^^{^k) wih exit from 
Ue(r]'') in some finite time ne(^/e,77*) and be attracted to a compact neighborhood B in 
some finite time le{^k)- Let b = d(C, B), then we have 

which is a contradiction to (Al). 

(2) If d{H'^''{^k)i K) < 6, then there exists some m^^ = m^{^i)^ 1 < z < m, such that 

which we go back to the first case, therefore, there is a contradiction to (Al). 

(3) If H^''{^k) ^ ^e(C*)5 then by Hartman-Grobman-Cushing theorem, wih exit from 
(7e(^*) in some finite time ^^(^/e^C*)? i-^-? 

From (Al), we have 

which we go back to either the first case or the second case, therefore, there is a 
contradiction to (Al). 

Based on the arguments above, we can conclude that the statement of Lemma A.l is true. □ 

Lemma A. 2. [Uniform weak persistence] If the system (3)-(4) satisfies Condition C1-C3. Then 
there exists some 5 > {), such that for any initial condition — (xo^yo) G X, the system has 

limsupx^ > 6. 

Proof. Since the system satisfies Condition C1-C3, then there exists a compact neighborhood 
W C Sy n S of the stable periodic-2 orbit My = {//i,'/72} attracting all points {0,y) G S from 
Theorem 4.3 (Elaydi and Sacker 2004). Condition C3 implies that My is transversal unstable, 
i.e., its external Lyapunov exponent is greater than 1. 
Define P{^) = x where ^ = (x, y) G X and 



Then r(t, •) is lower semicontinuous. For /i > 0, t > set 

a(h,t) = G X : r(t,0 > 1 + Z^}. 



Then a{h, t) is an open set from the semicontinuity and it has property that a{hi^t) C a{h2^ t) if 
hi > h2- Since W is compact, then there exists a /i > 0, and a finite increasing positive integers 
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ki^i = l...A^, such that 

N 
i=l 

where 

B = f^pi^ = e 1 : di^My) < e} 

is a compact neighborhood of W in X. We want to show that for any point ^ = (x, G S \ VF, 
its semi-orbit 7^(0 eventually exits from B with Xn > e for some positive integer n. If this is not 
true, then there exists some point ^ = (x,y) ^ B\W such that H'^{^) G S \ for all n > 0. 
Since any point in B belongs to some a(h^ ki)^l < i <h. This implies that there is a sequence of 
integer ni ^ oo with ni — E {fci, k^} for each i such that 

P {H^'^iO) > P (H'^'-'iO) [1 + /i] > ^[1 + hY ^oo since x > and ^ oc 

which is a contradiction to the fact that all points are attracted to the compact set D^. Thus, 
for any point ^ = G 5 \ VK, its semi-orbit 7^(0 eventually exits from B with > e for 

some positive integer n. Combined with Lemma A.l, we can conclude that for any point ^ G X 
that is close enough to C, it will enter the compact neighborhood B of My and exit from B in 
some finite time. Therefore, there exists some e > 0, such that for any initial condition ^ G X, 
the system has 

limsupx^ > e. 

□ 
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